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ABSTRACT 

Single white dwarfs with photospheric metal absorption and no dusty disks 
are believed to be actively accreting metals from a circumstellar disk of gas caused 
by the disruption and sputtering of asteroids perturbed by post-main sequence 
planetary systems. We report, for the first time, the detection of circumstellar 
Ca II gas in absorption around the metal-polluted white dwarf WD 1124-293 at 
>7 -Rwd an d <32000 AU, with a probable distance of ~54 Rwd- This detec- 
tion is based on several epochs of high-resolution optical spectroscopy around 
the Ca II H and K lines (A=3968A, 3933A) with the MIKE spectrograph on 
the Magellan/Clay Telescope. We confirm the circumstellar nature of the gas 
by observing stars with small angular separations to WD 1124-293 and larger 
distances from Earth-none of the reference stars show absorption consistent with 
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the presence of appreciable local interstellar medium at the distance of the white 
dwarf. By combining our observations over four years with previous spectra of 
WD 1124-293, we have measured the equivalent width of the two photospheric 
Ca lines over a period of 11 years. We see < 5-7% epoch-to-epoch variation in 
equivalent widths over this time period, and no evidence for long term trends in 
the strength of the Ca II lines. Since it is likely that the circumstellar gas arises 
from a disk edge-on to our line of sight, we also place limits to short period tran- 
siting substellar and planetary companions with R > R e using the WASP survey. 
The presence of gas in orbit around WD 1124-293 implies that most metal rich 
WDs could harbor relic planetary systems. Since roughly 25-30% of WDs show 
metal line absorption, the dynamical process for delivering small bodies close to 
a WD must be robust. 

Subject headings: circumstellar matter-planetary systems-white dwarfs 



Introduction 



In the absence of externally accreted material, radiative levitation, or convective dredge- 
up, white dwarfs (WDs) should have pure hydrogen or pure helium photospheres. All other 
elements heavier than helium should settle out of view on timescales much shorter than 
the cooling age of WDs. Hydrogen WDs (type DA) have photospheric settling times that 
can range from days to thousands of years, while helium WDs (type DB) have settling 
times as long as 1 Myr. Due to their high gravity and thin photospheres, WDs are some 
of the most sensitive probes to low levels of external accretion. A typical 0.6M Q hydrogen 
WD with an effective temperature of 10000 K can show observable Ca II lines from the 
accretion of ~10 6 g/s of solar abundanc e material, correspond i ng to the accretion of a small 



aster oid's worth of material per year (IPaquette et al.l Il986l ; iKoester et al.l 120051 ; iKoester 



20091 ). WDs are sens i tive probes o f the accretion from the surround ing interstellar medium 



(ISM) flDupuis et al 



Zuckerman et al. 



dust disks (e.g. 



1992 



1993a 



bl) , the winds from companions (jSion fc Starrfieldl 11984 



20031: lDebesll2006h. an d the accretion of asteroidal material from orbiting 



2010; Debes et al. 



2011a 



Zuckerman et al. 2007: Klein et al. 2010; Vennes et al. 2010; Dufour et al. 



b). All three of these phenomena can explain the presence of metal 



lines in WDs, but a significant fraction of metal polluted WDs (type DAZ for hydrogen 
atmosphere WDs and type DBZ or DZ for helium atmosphere WDs) are apparently single, 
do not po ssess any dust disks, and do not have sufficien t surrounding ISM to explain their 
accretion (lAannestad et al.lll993l ; iKilic fc Redfieldl 120071 ). 



ISM accretion was the first logical explanation for metal polluted WDs (jWesemael 



1979 
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Wesemael h Truranlll982l ; iDupuis et al.lll992l ll993aU bl). However, studies of polluted WDs 
show no strong correlation between v elocity and metal abundance as would be expected 



for Bondi-Hoyle accret ion of the ISM (lAannestad et al.l Il993l ; IZuckerman et al.l 12003k 12010 
Kilic &: Redfieldl 120071 ). N or is there any correla tion between metal abundance and con- 



straints on the local ISM (IKilic fc Redfieldl 120071 ). Most observed DAZs and DZs reside 
well within the low density Local Bubble, and are not physically located near known dense 
molecular clouds. The lack of appreciable hydrogen in DZs also cannot easily be explained 
by the accretion of solar abundance material. 

The presence of polluted DAZs that are members of WD+red dwarf binaries raises the 
possibility that previously unseen companio ns could be the cause of many polluted WDs 
through the accretion of a companion wind (IZuckerman et al.ll2003l ). While it is physicall y 
plausible that a close (a <1 AU) companion can pollute a WD with its wind ( IDebes! 120061 ). 
most WDs that have been observed with Spitz e r are constrained to have no close ste llar or 
brown dwarf companions (jMullally et al.l 120081 ; iDebes et al.l 120071 ; iFarihi et al.l 120091 ) . The 
most massive objects that could orbit these WDs would have masses of 5-15 Mj up . 

Dusty WDs all show pollution due to the direct accretion of metal rich grains onto the 
WD surface. However, only ~30% of metal po lluted WDs show infrared excesses due to 



dusty disks (IKilic et al.l 120081 ; IFarihi et al.l |2010| ) . Some disks may be close to edge-on to 



our line of sight and hidden, but inclination effects alone cannot explain the relative number 
of dusty and non-dusty WDs. Tidally disrupted planetesimals perturbed by planets are 



believ ed to be the origin of dusty WDs ( IDebes fc Sigurdssonl 12002 



Jura 2003; Bonsor et al. 



20111 ). Similarly, apparently single WDs with metal pollution may be caused by hitherto 
unseen disks that are either dusty, optically thin, and undetectable in the infrared, or entirely 
gaseous. Gaseous disks with some mass in gas of M disk that arise from the sputtering of debris 
from tidally disrupted planet esimals hav e been postulated to explain the relative number of 
dusty and non-dusty WDs (I Jural 120081 ). In this scenario, only planetesimals larger than 
~2 Mdisk create an optically thick dusty disk, while smaller objects are completely ablated 
by the relative velocity between the debris and the existing gaseous disk, further populating 
the disk with gaseous material. 

Circumstellar gas emission lines hav e been observed toward several hotter WDs that also 
possess dust disks (IGansicke et al.ll2006l . 120071 l2008t iMelis et al.ll2010l ). demonstrating that 
significant gas can be generated from either the sublimation, collision, or photoionization of 
dust grains. Similarly, a ha ndful of hotter WDs that have been observed in the UV show cir 



cumstellar gas absorp tion (IHolberg et al.l 1 1998c iBannister et al.l 120031 ; lLallement et al.l 12011 



Dickinson et al.l 120121 ). These hotter WDs may be evaporating circumstellar dust disks or 
local ISM, but to date none of these systems has been shown to possess a dust disk-their 
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origins are less clear. As WDs become hotter, radiative levitation of metals becomes a more 
important process and can explain t he presence of metals in a WD atmosphere for effective 



temperatures larger than ~25000 K (jChayer fc Dupuisll201Cll ). 



In this work we present the detection at 8-a of weak circumstellar Ca II absorption in 
orbit around the cooler DAZ WD 1124-293. Based on observations of other sight lines close 
to WD 1124-293 and the relative velocity of the circumstellar gas, we confine the spatial 
extent of the absorption to 7 R WD < r < 32000 AU from WD 1124-293. We argue that it 
must arise from a circumstellar gaseous disk in orbit around the WD at a probable radius 
of ~54 Rwd, well within the WD's tidal disruption radius. In §2] we describe the MIKE 
spectroscopic observations of WD 1124-293 and present the results of those observations in 
§21 We also search WASP data for planetary transits around WD 1124-293 in §4J We search 
for variability in the accretion of material on the WD in §5] and present our conclusions in 



Observations 



We have taken optical spe ctra of WD 1124-293 o ver 16 separate days spanning 4 years 
with the MIKE spectrograph ( IBernstein et al.l 120031 ) on the Clay Telescope at Las Cam- 
panas O bservatory between 2 007 and 2011 (See Tabled]) as part of a large DAZ monitoring 



project (IDebes fc Ki 



1C 



2010) . WD 1 124-2 93 was first discovered in the Edinburgh-Cape 



Survey and listed in iMcCook fc Sionl (119991) . Based on spectral m odels of WDs, WD 1124- 
293 has M=0.66 M and T eff = 9400 K koester fc Wilkenl lioo^i . with a spectroscopically 
determined distance of 33.6 pc ( iPauli et al.ll2006l ). Photospheric metal lines of Ca II were 
detected with mod erate resolution spectroscopy from Keck H IRES (R~30000) and VLT 



UVES (R-18500) f lZuckerman et al. 



2003t iKoester et al.l 120051 ). which showed moderately 



strong Ca II H and K lines (A=3968. 469, 3933. 664) with a heliocentrically corrected redshift 
of 29.5 km s^ 1 (29±1 km s" 1 ) for Keck (VLT) in December 1998 and April 1999 (April and 
May 2000). The calcium cannot be explained by accretion of the local ISM based on mea- 
surements of the local ISM density ( jKilic fc Redfieldl 120071 ) . The equivalent width (EW) of 
the Ca II line over the four separate epochs with VLT and Keck averaged f=al06 mA with no 
significant variation over this period, implying a steady state accretion rate of 4xl0 15 g/yr 
( IKoester fc Wilkenl 120061 ) . assuming 1% metal rich material relative to solar abundance. In 
January and June of 2010 we observed four bright F and A stars located near WD 1124-293 
on the sky in order to probe the surrounding ISM at the distance of WD 1124-293. 

For our MIKE spectra, we used a 0'/7x5" slit which corresponded to a resolution of 
~40,000 at the Ca II K line. Th-Ar comparison lamp spectra were taken near in time to 
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each spectrum of WD 1124-293. In addition, we obtained the previously reduce d spectra of 
this WD from the VLT and Keck faoester et al.ll2005l : IZuckerman et al.ll2003h . The VLT 
spectra were kindly provided by D. Koester, while the Keck data was publicly available with 
the journal article it was published in. 



Kelson et al. 


(2000) and 


Kelson 


(2003) 



We initially extracted the spectra using the full pipeline, but discovered that the pipeline 
was subtracting the background at a level of precision insufficient for high quality equiv- 
alent width measurements. To rectify this, we performed our own extraction, background 
subtraction, and wavelength solution for 8-10A windows centered on the Ca H and K lines 
for all of our MIKE spectra. The spectra were also corrected for heliocentric velocity shifts 
between observations. Since the Ca II K line falls on two orders, we combined each order by 
fitting the continuum with a polynomial and averaging the orders for each night to increase 
signal-to- noise (S/N) to the values reported in Table [TJ All the fitted spectra from MIKE 
were then averaged with a S/N weighting to produce a final spectrum of the region around 
each line. 

In order to check for variability in our spectra between 2007-2011 as well as with the 
previous observations, we calculated the EW of the two calcium lines from each night's 
spectrum. To calculate the EW, we chose a window around the line equivalent to ±3 times 
the standard deviation of the line as determined by a Gaussian fit, which also returned 
a measure of the radial velocity for the WD on each night. We chose several different 
polynomial fits to the continuum and added any systematic uncertainty in the continuum fit 
to our uncertainties in the EW calculation. This approach has the advantage of calculating a 
more accurate EW despite differences in resolution between different epochs, since the total 
integrated line flux is conserved from observation to observation. Table [2] details the observed 
photospheric line radial velocities and equivalent widths. From the individual measures of the 
K line with MIKE, we obtain a median radial velocity for WD 1124-293 of 29.0±0.5 km s" 1 . 



3. Results 

The final, S/N weighted spectrum of WD 1124-293 at the Ca H and K lines is shown in 
Figure [TJ Just to the blue of the main photospheric lines we detect weak lines indicative of 
another Ca II component, with line depths of 3.4- and 8.5-cr below the continuum for the H 
and K line respectively. As seen in Table |2l the lines at both wavelengths are separated by 
^30 km s _1 to the blue of the photospheric Ca II. Weak lines in the spectrum could be due 
to additional photospheric species, circumstellar Ca II gas, intervening ISM Ca II gas, or an 
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artifact from the spectrograph or reduction process. In this section we argue that the lines 
arise from circumstellar Ca II gas in orbit around the WD. 



3.1. Ruling out artifacts and other atomic species 

There are three separate conditions that rule out the possibility that the lines could 
arise from the data reduction process or the MIKE spectrograph. First, we have verified that 
the line near the photospheric Ca II K line was detected significantly on a S/N combined 
spectrum of each of the two echelle orders that the line fell on. The position of the line on the 
two orders is separated by ~1100 pixels, ruling out possible artifacts from the instrument. 
For both orders, the line is detected at the proper position and at 4- and 6-cr below the 
continuum. Secondly, the line is detected significantly in the final combined spectra when 
using the full Kelson MIKE pipeline, which used different methodology to perform order 
extraction and background subtraction. Finally, when we combine the data into subsets, we 
recover the line as well. Figure [2] shows three separate epochs: March 2008, May 2008, and 
a combination of our June 2010 and June 2011 data. In particular, the heliocentric velocity 
correction in March 2008 and our June 2010+2011 data differed by 23 km s _1 , ruling out 
any terrestrial origin for the line, since the epoch-to-epoch centers match to within 0.01A. 

These lines could be weak photospheric features. We searched atom ic line lists, as we ll 



as metal lines other than calcium reported for other polluted WDs (e.g. iKlein et al.ll2010l ). 
and found no atomic lines that could plausibly be matched to the observed spectrum of 
WD 1124-293. Furthermore, the line depths and relative equivalent widths are consistent, 
within the uncertainties, to the expected line strength ratio between the Ca II H and K 
lines. Based on these arguments, the weak features must arise from intervening ISM or 
circumstellar Ca II gas. 



3.2. Ruling out Local Interstellar Absorption 



While ISM absorption could explain the presence of the weak features we observe, we 
rule out an ISM contribution. The Sun resides in a low-density ISM structure known as 
the Local Bubble. WD 1124- 293 is at a spect roscopically determined distance of 33.6 pc, 
well inside the Local Bubble (jPauli et al.ll2006l ). The local interstellar medium is relatively 
devoi d of cold neutral gas, until the accumulation of dense material at ~80 pc (IWelsh et al. 



2010). However, warm, partially ionized clouds have been detected throu gh weak Ca II ISM 



absorption features in the spectra of stars within 50 pc of the Sun (see iRedfield fc Linsky 
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20081 ; IWelsh et al.ll2010l . for a review of the properties of the local interstellar medium). 



If the weak Ca II absorption feature in WD 1124-293 is due to ISM accretion, other 
stars in close proximity on the sky should also show a similar absorption feature. In order to 
test this possibility, we observed several bright stars at angular separations < 1.6° from WD 
1124-293 to search for absorption features at the same velocity. This technique, as well as 
searching for line variability, has been used in the past to co nfirm circumste llar absorption 



around HD 32297, a young star with an edge-on debris disk (IRedfieldl 120071 ). We observed 
four stars with constraints on their distance from Hipparcos and with the same setup as WD 
1124-293: HIP 55864 (F5V, r=0.27°) HIP 55901 (A0V, r=0.71°), HIP 55731 (A9V, r=1.10°), 
and HIP 55968 (A3V, r=1.58°) (See Tabled). Another star, HIP 56280A (F8V, r=1.19°), 
had a spectrum that was kindly provided to us by S. Desidera. Figure [3] shows how each 
ISM standard is arranged on the sky relative to WD 1124-293. 



The distances for each star are given in the latest reduction of Hipparcos data (Ivan Leeuwen 



20071 ). HIP 56280A and HIP 55864 have secure distances (parallax measurements with er- 
rors of <20%), while HIP 55731 and HIP 55901 have marginally detected parallaxes (~3-er). 
HIP 55968's parallax is consistent with zero. For the marginal detections, we take as the 
stellar distance the 3-a upper limit to the parallax, which is 6.62 mas (151 pc) and 5.67 mas 
(176 pc) for HIP 55731 and HIP 55901 respectively. HIP 55968's optical and NIR magnitudes 
and spectral type are consistent with this star lying between HIP 55731 and HIP 55901, so 
we take the lower limit to its distance to be 151 pc. 

Each star had broad Ca II K photospheric absorption which we removed by apply- 
ing a polynomial fit to the photospheric line, an example of which is shown in Figure HI 
For HIP 56280A, its narrow line core was fitted with two polynomials joined at the line 
core center of 3933. 704A. The resulting continuum fitted spectra were inspected for narrow 
absorption components. HIP 56280A (to lower significance) and HIP 55864 show no appre- 
ciable absorption features comparable in strength to the WD 1124-293 line, while HIP 55731 
shows a weak (presumably ISM) component far from the radial velocity or our observed 
feature. HIP 55901 and HIP 55968, both constrained by Hipparcos to be more distant than 
HIP 55864, each showed absorption from at least three ISM components. We attempted 
to simultaneously fit these three components for the two stars and present our best fits in 
Figure [5j which shows WD 1124-293's Ca II K line region compared with our ISM standards. 
Table [3] shows the resulting velocity components for these fits, consistent with at least three 
distinct complexes of ISM gas that slowly vary on size scales of ~l-2 pc, one of which is at 
a radial velocity within a few km s _1 of the circumstellar feature. Despite the similar radial 
velocity of the two features to what is observed in WD 1124-293, the non-detection of any 
similar lines for HIP 56280, HIP 55864, and HIP 55731 constrains the location of these ISM 
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components to be at a distance >151 pc. The non-detection of the faint candidate circum- 
stellar component for HIP 55864 constrains the extent of the Ca II gas near WD 1124-293 
to < 0.27° at d=33.6 pc, or roughly 32000 AU (0.16 pc). 



Welsh et al.l ( 120101 ) find a slowly increasing Ca II equivalent width with increasing dis- 



tance, but depending on whether local gas clouds are encountered or not, they find both low 
and high va l ues of Ca II volume density within 30 pc. Out of the 1857 stars analyzed by 
Welsh et al.l (120101 ) . 50 are within 40 pc of the Sun. Only one of these 50 stars, HD 159561 
(I = 36°, b = +23°), shows an ISM absorption feature with a Ca II equivalent width stronger 
than what is observed for the weak feature in WD1124-293 (/ = 282°, b = +30°). The ISM 
complexes observed for our more distant stars are consistent with a rel atively dense clou d 
of material that lies >90 pc along the Galactic sightline to these stars (jWelsh et al.ll2010l ). 
and thus cannot account for the weak Ca II feature we measure for WD 1124-293. Hence, 
it is unlikely that this feature is due to an intervening warm ISM cloudlet. The only sce- 
nario that the data allows is if an extremely narrow warm ISM column near WD 1124-293 
(width ~0.06 pc and height ~1.3 pc) with a velocity gradient intersected the sightlines to 
WD 1124-293, HIP 55901, and HIP 55968, but did not intersect the sightline of HIP 55864. 



Such a pathological structure is not supported by studies of the nearby ISM. lRedfield fc Linsky 



( I200ll ) studied the small scale (0.05-1.2 pc) structure of the local interstellar cloud and find 
that the Mg II column densities do not vary by more than a factor of two for < 0.6 pc 
scales. The implied Ca II column densities by the much stronger ISM absorption observed 
towards HIP 55901 and HIP 55968 are an order of magnitude larger than for WD 1124-293, 
and would have been easily detected in the ISM standards that showed no component at the 
velocity of our candidate circumstellar feature. Hence, the observations of our ISM standards 
limits the amount of interstellar Ca II between WD 1124-293 and us at the velocity of the 
Ca II K feature to EW < 6 mA. 

The source of metals in high surface gravity WDs has been a puzzle for decades (see 
Dupuis et al.lll993al ) . The search for a corr elation between the loca l I SM clouds and the loca- 
tions of metal-rich WDs was inconclusive (jZuckerman et al.ll2003l ). iKilic &: Redfieldl (120071 ) 
used the observed ISM column densities toward stars in close proximity to the known DAZ 
WDs to demonstrate that there is no correlation between the accretion density required to 
supply metals observed in DAZs with the densities observed in their interstellar environment, 
including WD 1124-293. Therefore, the DAZ WD population as a whole argues against ISM 
accretion being the dominant contributor to the metals in WDs. 

In conclusion, it is more likely that the weak Ca II feature is due to circumstellar 
gas in orbit around WD1124-293 because (1) no similar absorption features are detected 
in neighboring sightlines as close as 0.16 pc and the local interstellar medium is relatively 



-9 - 



homogeneous at such scales (IRedfield fc Linskyl l200ll ). (2) the strength of the absorption 
feature is uncommon for stars within 40-50 pc of the Sun, (3) the absorption feature matches 
the stellar radial velocity, and (4) there i s no correlation betwe en the metal abundances of 
known DAZ stars and the ISM densities (jKilic fe Redfieldl 120071 ) . 



3.3. Limits to the inner extent of circumstellar gas 



If the line arises from gas in orbit around the WD, the relative velocity between the pho- 
tospheric line and the circumstellar feature can be used to determine how deep the gas lies 
in the gravitational well of the WD. The photospheric Ca II line is offset from the true sys- 
temic velocity of the WD due to the gravitational redshift of WD 1124-293. Using the latest 
values of WD 1124-293's log g=8.096 and its T cff =9420±150 K (D. Koester, private com- 
munication), we have calculated masses from cooling; models of WD jj[ ( iBergeron et al.lll995l 



Holberg fc Bergeron 



radius relation of 



2006 



Fontaine et al. 



Tremblav et al.l 1201 ll ; iBergeron et al.l 1201 ll ) which use the mass- 



fbooiT I" to determine the gravitational redshift of WD 1124- 
293. The uncertainty in the gravity of the WD dominates the uncertainty in the expected 
gravitational redshift, but is also no t well known. To esti mate this poss i ble u ncertainty we 
compare our value to that given in Koester et al.l ( 120011 ) and iKoesterl ( 120091 ). which vary 
between 8.04 to 8.099. We therefore estimate the systematic uncertainty in the WD gravity 
to be ps0.06. Using this value and uncertainty and calculating a gravitational redshift (v : 



grav J i 



we obtain t> grav =34.9±3.7 km s 1 , implying a systemic velocity 7=-5.9±3.7 km s 1 



Previous calculations of WD 1124-293's 7 (based on its inferred mass and r adius from 
synthetic sp ectral modeling) find -2.7±4.3 km s -1 , consistent with our calculated 7 (jPauli et al 
2003L 12006( 1 . The primary difference between our result and previous calculations of 7 come 
from a higher WD gravity (The earlier values of log g being 0.05 dex lower), highlight- 
ing the dependence of this calculation on a detailed knowledge of a WD's properties. The 
measured relative velocity between the photospheric line and the circumstellar line (v gas ) 
from our S/N weighted spectrum is -29.9±0.8 km s _1 . The orbital radius of the gas is then 
R g =v gra _ v / (f gas + fgrav) Rwd, implying a minimum radius to the gas of 7^ Rwd- Given the 
uncertainties, the faint line could correspond to the systemic velocity and reside further from 
the WD, thus directly probing the gravitational redshift of the WD. This minimum radius 
then should be treated as a low er limit to t he gas radius. For optically thin silicate dust, the 
sublimation radius is given by (I Jural l2008t iKimura et al.l 120021 ; lHoardll201ll ): 



These models are available at 



http://www.astro.umontreal.ca/~bergeron/CoolingModels 
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V 0.001 l q> / V ^ub J 



R suh = 3.7 x 10+ u - cm. (1) 



From the cooling models used to determine WD 1124-293's gravitational redshift, we 
also determine a luminosity of almost exactly 0.001 L , leading to a sublimation radius of 
~44 -RwD; which is close to the inferred radius of the gas from kinematic arguments in §3.41 



3.4. Properties and Origin of the Circumstellar Gas 

Further characterization of the Ca II lines permits additional constraints on the origin 
and properties of the circumstellar gas around WD 1124-293. In particular, the different 
constraints point to a localized disk of gas a few tens of WD radii from WD 1124-293 and 
well within its tidal disruption radius. 

In order to determine the column density of the Ca II gas, we minimized a \ 2 metric 
for an optically thin gas column by varying the column density (log N) and FWHM of the 
line, using the proper oscillator strengths and assuming a Gaussian line profile. From the 
minimization we find a best fitting line model with (98% confidence) log iV=llio2 cm -2 and 
FWHM=0.2±0.1 A with a x^=l-08. The best fit line, compared to the observed circumstellar 
line is shown in Figure O The model overpredicts the expected line strength of the Ca II H 
line, but given the low S/N of the line it is hard to determine if this difference is significant. 
Since this is the ionized state of calcium, it is impossible to tell how much total gas resides 
in the system without a) other gas species or b) a measurement of Ca I to determine the 
temperature and electron density of the gas. Further characterization of the gas will be the 
focus of future work, but is beyond the scope of this paper. 

We can also constrain the location of the gas by the width of the weak Ca II line-an 
upper limit to the orbital velocity of the gas can be determined from the line FWHM, and 
thus a characteristic radius for the disk. We assume that the gas is in a circular orbit and that 
the line FWHM corresponds to the range of orbital radial velocities within the circumstellar 
gaseous disk that intersects the angular disk of the WD on the sky: 

?"disk ~ (GMwd) 1 

From our gravitational redshift calculations, -Rwd ~8.4x10 8 cm, and the FWHM of the 
line is 16 km s -1 , giving a radius for the disk of 4.5 xlO 10 cm or «54 .Rwd- This is consistent 



2R 



WD 



FWHMicm s~ v 



(2) 
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with the gas being close to the WD and within the tidal disruption radius of WD 1124-293, 
as well as being from the sublimation of optically thin dust. 

The constraints on the location of the gas and the presence of Ca II combined with 
the lack of any appreciable ISM all point to a similar origin for WD 1124-293 as that 
posited for WDs with near-IR excesses due to dusty disks. These observations fit in with 



namically perturbed by a planetary system ( 


Debes & Sieurdsson 


2002: 


Bonsor et al. 


2011: 


Bonsor & Wvatt 


2012: 


Debes et al. 


2012 


) and then tidally disrupted 


(Jura 


2003 


) . How- 



ever, a size distribution of planetesimals will be disrupted over time, and for smaller bodies, 
sputtering renders most tidally disrupted dust into a gaseo us phase quite quickly, formin g 



circumstellar gas disks instead of circumstellar dust disks (IJural 120081 ; iFarihi et al.l 120121 ) . 
Circumstellar gas absorption roughly at the same radius as that observed for dusty disks 
around WDs is also suggestive of a tidally disrupted asteroidal origin for our observations of 
WD 1124-293. 



4. Limits to Transiting Companions to WD 1124-293 

Given the possible edge-on inclination of this system due to the presence of circumstellar 
gas, the probability of any detectable planetary companion transiting WD 1124-293 might be 
higher than a randomly inclined system. To that end we have gathered all data of WD 1124- 
293 from the Wide Angle Search for Planets (WASP) transit database. From previous 
observations, WD 1124-293 does not possess close-in brown dwarf or stellar companions. In 
their comprehensive Spit zer search for IR excesses due to companions around a sample of 



WDs, IFarihi et al.l ( 120081 ) placed a definitive upper limit of 12 M Jup to the presence of any 



unresolved companions in orbit around WD 1124-293. 

The WASP North and South telescopes are located in La Palma (ING - Canary Islands) 
and Sutherland (SAAO - South Africa), respectively. Each telescope consists of 8 Canon 
200 mm //1.8 focal lenses coupled to e2v 2048x2048 p ixel CCDs, yielding a field-of-view 



of 7.8x7.8 square degrees with a pixel scale of 13.7" (jPollacco et al.l 120061 ) . The WASP 
observing strategy is such that each field is observed with a typical cadence of the order of 
8-10 minutes and typical exposures of 1 minute (30 sec exposure plus over heads). WASP 
provides good quality photometry with a precision exceeding 1% per observation in the 
approximate magnitude range 9 < V < 12. 

WD 1124-293 was observed by the WASP-South telescope in the 2007 and 2008 seasons 
covering the interval 2007 January 04 to 2007 June 05, and 2008 January 05 to 2008 May 
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28, respectively. The 10618 pipeline-processed photometric measurements were de-trended 
using the Tamuz algorithm (ITamuz et al.ll2005l ) to account for linearly-correlated systematic 
noise in the data. We show the WASP light-curve in Figure [7] (bottom-panel). 

We u sed the modified implementation of the Box-Least Square (BLS) algorithm de- 
scribed in iFaedi et al.l ( 120111 ) to search for transits and eclipses of sub-stellar and planetary 
companions in close orbits around the white dwarf WD1124-293. We searched a parame- 
ter space defined by orbital periods ranging from 2 hours to 15 days, and companion sizes 
ranging from Moon-size to Jupiter-size. No transiting planet has been identified. Figure 
(top-panel), shows the BLS power-spectru m for the sta r WD1 124-293. The dashed-line indi- 
cates the detection threshold as defined in IFaedi et al.l (1201 if ) for 10% noise levels. We have 
phase-folded the WASP light-curve to investigate the possible presence of an astrophysical 
signal at the periods detected above the threshold. However, none showed real variation. 
Finally, we have used the modified Lomb-Scargle periodogram ( jScargld 119821 ) to search for 
generic variability. We find that a real power is assigned to a peak at P=27.42 days which 
is most probably related to the Moon cycle. We used these null results together with the 
results of our simulations to estimate an upper-limit to the frequency of close companions 
to WD 1124-293. 

To assess the chances of detecting eclipses and transits of substellar and planetary 
companions to WD1124-293, we performed an extensive set of Monte Carlo simulations. 
The approach we adopted was to create realistic synthetic light-curves containing eclipse 
and transit signatures of the expected depth (5 tr ) and duration for a range of companion 
sizes and orbital periods, then to attempt to detect these signatures using the modified BLS 
mentioned above. By noting the rate at which the BLS search recovered the transit (fdet) at 
the correct period (or an integer multiple or fraction) we were able to estimate the detection 
limits of such systems in an automated manner. 

Details of the synthetic transit signal g enerated, the corre sponding transit probability, 
depth and duration are described in detail in lFaedi et al.l (120111 ) . Here we present the results 
of simulations covering the orbital period-planet radius parameter space defined by seven 
trial periods spaced approximately logarithmically (P = 0.08, 0.22, 0.87, 1.56, 3.57, 8.30 and 
14.72 days), and five planet radii R p = 10.0, 4.0, 1.0, 0.55, and 0.27 R e . 

Figure [8] shows the detection rate of the simulated transit signals injected in the WASP 
light curve of WD1 124-293. We regard as a match any trial in which the most significant 
detected period is within 1% of being an integer fraction or multiple from 1/5 x to 5x the 
injected transit signal. From our result on the simulations conducted using the WASP light 
curve of W D1124-293 (Figur e El see also results on simulations conducted on WASP data 
detailed in lFaedi et al.l 1201 lh . we would have been able to detect the transit signal due to 
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a companion of radius lR ffi and period between 2 and 16h with a confidence of 75%. For 
objects of radius similar to that of the Earth but slightly bigger (e.g. up to 3R ffi , or so-called 
SuperEarths) we would have been able to detect the eclipse signal for periods between 2h 
and 12h with a confidence of 95% and to about Id with 75% confidence. Eclipsing gas 
giant planets and substellar objects would have been detectable at orbital periods between 
2h and 1.7d with 95% confidence, and with orbital periods of 2d with 75% confidence. On 
the other end, no stringent constraint can be placed on objects with radius < 1R® even at 
short periods. The same is true even f or large bod ies such as Jupiter with orbital periods 



longer 2 days. From Figures 1 and 2 of lAgoll (120111 ) and using the white dwarf temperature 
of 9420 K and mass of 0.66 M Q we find that for WD1124-293 the white dwarf habitable 
zone (WDHZ) extends for a region comprised between 0.02 and 0.045 AU for which, in our 
simulations, we can only put weak constraints for Neptune/Uranus-size (~ 3 — 4R e ) and 
larger objects. However, eventually WD 1124-293 will cool further and when it cools to below 
7000-9000 K, planets with orbital separations of 5xl0~ 3 - 0.02 AU (P ~0.16-1.3 days) will be 



in a "continuous habitable zone" which lasts for >3 Gyr (lAgoll [20 111) . Our observations thus 
rule out a significant fraction of > Earth-sized planets in this system that might reside in 
the continuous habitable zone of WD 1123-293. 

We investigated the various factors which affect the efficiency of detection of these transit 
signals. When generating each synthetic light-curve we can easily assess a priori whether it 
will fail the tests requiring a minimum number of individual transits (> 5) and in-transit data 
points (> 5). We find that our requirements alone render transiting companions essentially 
undetectable at our longest trial periods (8.30 and 14.72 days); the transits are too short 
in duration and too infrequent to be adequately sampled. For companions around 1R® and 
larger however there is a good chance of detection out to periods of ~2 days. 



Limits to a variable Accretion Rate to WD 1124-293 



Given the stochastic origins to dust around WDs, it might be possible to observe variabil- 
i ty in WD accretion rates. Ca v ariability was claimed to be present in t he dusty DAZ G 29-38 
Jvon Hippel fc Thompson 2007 ). but disputed with other observations ( Debes fc Lopez-Morales 
20081 ). Discrepant abundances of Si between optical and UV obse rvations of some dusty 
WDs may suggest possible variable accretion (IGaensicke et al.ll2012l ). The multiple spectra 
of WD 1124-293 allow a careful look at the variability in the accretion rate of material onto 
the surface of the WD through EW measurements of its photospheric Ca II lines. WD 1124- 
293's metal settling timescale for calcium is 194 yr, much longer than the total time high 
quality spectra of the Ca line have been available, but short enough that small changes due 
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to settling may be observed. Similarly, if there were any short duration accretion events that 
briefly increased the accretion rate, we may have observed a sudden increase in the Ca EW. 

Figure |S] shows the EW as a function of Julian date (JD). Two measurements of the 
Ca II K line, from May 2009 and July 2010, differ significantly from the other measurements. 
Similarly, the Ca II H line measurements from 1998 and March 2007 also deviate significantly. 
The cause of these discrepancies most likely due to lower S/N and poor seeing, but given 
the good agreement between the other EW measurements and the lack of similar variability 
at the same time from the complementary Ca II line, we neglect these measurements in 
our analysis. To determine whether WD 1124-293's Ca II line strength was varying, we 
calculated a x 2 value for a constant value EW, determined from the median value of the 
observed EWs. With this measure, the Ca EW has a x 2 value of 23 (14) for the Ca II K (H) 
line, which corresponds to a probability th at the EW is constan t to 0.03 (0.3). Based on our 



full survey of other DAZs for variability ( iDebes fc Kilid 120101 ). we require the probability 
of constant EW to be less than 10~ 3 to be considered significantly variable. The standard 
deviation of our MIKE measurements, neglecting the anomalous measurements, is 4.5% for 
the K line and 7% for the H line. 

The roughly constant flow of material over ~11 years for this system places a firm 
lower limit on the amount of metal-rich material recently deposited, 4.6 xlO 16 g, once again 
assuming metal rich material at 1% the solar abundance values, and a to tal amount of 



mate rial over one settling time of 9x10 g assuming a settling time of 194 yr ( IKoester et al. 



20051 ). These masses correspond to asteroids with sizes of 1.5 and 9 km respectively, assuming 



a bulk density of 3 g/cm 3 . These radii are comparable to the sizes of a steroids th at should 



survive post-main sequence evolution at an orbital radius of a few AU ( 1 Jural 120081 ). 



6. Conclusions 

We have obtained moderate resolution optical spectroscopy of the DAZ WD 1124-293 
over several epochs and found a weak circumstellar absorption feature of Ca II consistent 
with a tenuous gaseous disk supplied by small rocky bodies. Steady state accretion of 
material onto the WD surface appears to be occurring despite the relatively longer inferred 
settling timescale for metals, consistent with a recent disruption of an asteroidal body that 
was relatively small. 

WD 1124-293 has also been searched for transiting companions given that any planetary 
system is nearly edge-on to our line of sight, though none are found for short periods. Our 
results do not rule out the possibility for planetary companions at wider separations. 



- 15 - 



Emission from gaseous disks have be en observed around metal enriched WDs with disks 



that have been actively accreting dust (IGansicke et al.l 12006 



2007, 



2008). The origin of 



that gas could arise fr om the viscous h eating of material (jWerner et al. 2009 ) or from the 
photoionization of gas ( IMelis et al.ll2010f ). However, these systems may be unusual in the fact 
that they may represent an early stage in the creation of a dusty disk or represent a massive 
example of a dusty disk in orbit around a WD. The accretion of many small planetesimals 
shoul d proceed with little observable dust, and a primarily gaseous circumstellar disk (IJura 



20081 ). These properties match WD 1124-293, representing the first example of a single metal 



enriched WD with a direct link to the tidal disruption of planetesimals and thus to a relic 
planetary system. The large fraction of single WDs that show metal pollution implies a 
large fraction of post-main sequence objects may house remnant planetary systems with a 
robust mechanism for delivering small objects in WD-grazing orbits. A large spectroscopic 
survey of apparently single, dust-free WDs with high S/N may reveal more systems similar 
to WD 1124-293. 



We wish to thank Detlev Koester for unceasing aid in discussing the fundamental pa- 
rameters of WD 1124-293, and Silvano Desidera for kindly providing valuable spectra of 
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ing of spectral absorption lines, which provided a key impetus for constraining the location 
of WD 1124-293's gaseous disk. Finally, we want to thank the anonymous referee, whose 
attention to detail greatly increased the quality of this paper. This work made significant 
use of the SIMBAD and VIZIER databases, operated at CDS, Strasbourg, France. 
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Table 1. Spectroscopic Observations 



UT Date (JD-5400000) 


Target 


Observing Setup 


S/N/pixel (@Ca II) 


Notes 


1998-12-12 (51159.629) 


WD 1124-293 


Keck/HiRES 


37 


R=30,000 


1999-4-19 (51288.328) 


WD 1124-293 


Keck/HiRES 


26 


R=30,000 


2000-4-23 (51657.663) 


WD 1124-293 


VLT/UVES 


36 


R=18500 


2000-5-19 (51681.588) 


WD 1124-293 


VLT/UVES 


14 


R=18500 


2007-3-31 (54189.682) 


WD 1124-293 


Clay/MIKE 


24 




2008-3-21 (54546.605) 


WD 1124-293 


Clay/MIKE 


50 




2008-3-22 (54547.656) 


WD 1124-293 


Clay/MIKE 


50 




2008-5-13 (54599.536) 


WD 1124-293 


Clay/MIKE 


60 




2008-5-14 (54600.551) 


WD 1124-293 


Clay/MIKE 


49 




2008-5-15 (54601.477) 


WD 1124-293 


Clay/MIKE 


52 




2008-5-16 (54602.502) 


WD 1124-293 


Clay/MIKE 


56 




2009-4-16 (54937.686) 


WD 1124-293 


Clay/MIKE 


39 




2009-5-18 (54969.548) 


WD 1124-293 


Clay/MIKE 


38 




2010-1-1 (55197.727) 


HIP 55731 


Clay/MIKE 


36 


ISM Standard 


2010-1-1 (55197.733) 


HIP 55864 


Clay/MIKE 


40 


ISM Standard 


2010-1-1 (55197.738) 


HIP 55901 


Clay/MIKE 


100 


ISM Standard 


2010-1-1 (55197.743) 


HIP 55968 


Clay/MIKE 


47 


ISM Standard 


2010-6-17 (55365.467) 


HIP 55731 


Clay/MIKE 


78 


ISM Standard 


2010-6-17 (55365.472) 


HIP 55901 


Clay/MIKE 


162 


ISM Standard 


2010-6-17 (55365.479) 


HIP 55864 


Clay/MIKE 


105 


ISM Standard 


2010-6-17 (55365.486) 


WD 1124-293 


Clay/MIKE 


73 




2010-6-17 (55365.575) 


HIP 55968 


Clay/MIKE 


56 


ISM Standard 


2010-7-07 (55384.500) 


WD 1124-293 


Clay/MIKE 


33 




2011-6-09 (55385.493) 


WD 1124-293 


Clay/MIKE 


70 
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Table 2. Properties of Photospheric and Circumstellar Ca II lines in WD 1124-293 



Line type 


JD-5400000 


Ca II K EQW 


Ca II H EQW 


r*\ TT TV" 

Ca II K v r 


TT TT 

Ca II H v r 






(mA) 


(mA) 


(km s 


(km s J 




51159 629 


104± 3 


39± 4 


30.3 


31.3 




51288 328 


98± 5 


57± 6 


28.7 


27.5 




51657 663 


111± 6 


58±12 


29.1 


28.0 




51681 588 


109±17 




29.8 






54189 682 


106± 4 


34± 4 


29.7 


34.5 




04040. DUO 


1 1 «4_ 


0z± o 




zy.o 




54547.658 


113± 2 


54± 3 


29.0 


29.4 




54599.536 


113± 2 


53± 2 


28.3 


28.3 




54600.551 


114± 3 


49± 3 


28.6 


28.8 




54601.477 


113± 2 


55± 2 


29.0 


30.1 




54602.502 


118± 2 


55± 2 


29.2 


29.4 




54937.686 


113± 3 


53± 3 


28.8 


27.7 




54969.548 


127± 3 


43± 3 


29.1 


27.8 




55364.486 


115± 2 


52± 1 


28.7 


29.5 




55385.493 


95± 3 


59± 4 


28.4 


29.7 




55721.598 


111± 2 


54± 2 


29.7 


29.4 


Circumstellar 


March 2008 


11±2 




-0.5±0.9 






May 2008 


12=1=1 




-2±1 






June 2010+ June 2011 


10±1 




-1±1 






Total 


11±1 


2±0.5 


-0.9±0.5 


-0.4±0.5 
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Fig. 1. — Detection of weak circumstellar Ca II features in the spectrum of WD 1124-293. 
The spectrum is from a S/N weighted combination of observations spanning four years, and 
has a S/N of 180. The photospheric line is located at ^3934A and the circumstellar feature 
is blue shifted from the photospheric line. 
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Fig. 2. — Three separate epochs which show a clear detection of the blueshifted weak Ca II 
feature. The line is recovered relative to the continuum with a significance of 4.7-, 4.9-, and 
4-ct for March 2008, May 2008, and June 2010+June 2011 respectively. 
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Fig. 3.— Location of ISM standards on the sky relative to WD 1124-293. HIP 55864 is the 
closest star to WD 1124-293 and is at a projected separation of 0.16 pc at the distance to 
the WD. 
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Fig. 4. — Comparison of polynomial fitting of the broad and deep photospheric Ca II lines of 
two ISM standards, HIP 55731 and HIP 55864. The histograms are the observed spectra in 
each order of the MIKE echelle spectrum while the solid lines are the fits. The solid vertical 
line marks the location of WD 1124-293's circumstellar feature. 
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Fig. 5. — Comparison of the location of the weak Ca II K feature in WD 1124-293 with our 
ISM standards, arranged by probable distance from Earth. The dark vertical line shows the 
location of the weak circumstellar feature of WD 1124-293, while the horizontal blue line 
shows the 3-cr detection threshold for standards that showed no absorption at the wavelength 
of WD 1124-293's component. More distant standards show ISM absorption, while standards 
closer than 124 pc show no appreciable ISM absorption, ruling out LISM as a source for 
WD 1124-293's blueshifted feature. 
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Fig. 6. — Comparison of the observed circumstellar Ca line (black histogram with red error 
bars) with the photospheric line divided out to a model Ca II cloud (solid dark line) with a 
column density of log iV=11.07!ta2- 
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Fig. 7. — Top panel: Box-Least Square power spectrum for the white dwarf WD1124-293. 
Lower panel: WASP light-curve for WD1124-293 covering the period January 2007 - May 
2008, with a total of 10618 photometric measurments. We used WASP V-magnitude, where 
V mas = —2.5 * log(f lux/ 10 6 ), expressed in //Vega. 
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Table 3. Properties of ISM clouds for ISM standards 



Name Cloud I Cloud 2 Cloud 3 Cloud 4 

v EW v EW v EW v EW 

(km s^ 1 ) (inA) (km s _1 ) (mA) (km s^ 1 ) (mA) (km s -1 ) (mA) 



HIP 56280 




<10 




<10 




<10 


<10 


HIP 55864 




< 7 




< 7 




< 7 


< 7 


HIP 55731 




< 6 




< 6 


8.7 


22±2 


< 6 


HIP 55901 


-10.3 


47±1 


-1.7 


34±1 


9.5 


17±1 


<2 


HIP 55968 


-6 .3 


72±2 


2.0 


43±2 


11.8 


56±2 


21.2 19±2 




Fig. 8. — Detection limits for planetary companions of sizes ranging from Jupiter-size to 
Moon-size in orbit around the WD1124-293. The numbers on the contours relate to the 
detection efficiency in percentage. 
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Table 4. Recovery rate of simulated transits of WD1124-293 



Size 


R p 


5tr 


P 


fdct 




ft© 


(%) 


(days) 


(%) 


Jupiter 


10.0 


100 


0.08 


100 








0.22 


100 








0.87 


100 








1.56 


91 








3.57 


8 








8.30 





Uranus 


4.0 


100 


0.08 


100 








0.22 


100 








0.87 


99 








1.56 


82 








3.57 


1 








8.30 





Earth 


1.0 


49 


0.08 










0.22 


99 








0.87 


61 








1.56 


1 








3.57 
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Fig. 9. — The EW of Ca II in WD 1124-293's photosphere as a function of time. The EWs are 
compared to the expectation of a constant value, and do not significantly vary. Neglecting 
spurious points, the largest epoch-to-epoch variations are less than 7% and 5% from the 
median EW value for the Ca H (lower points) and K (upper points) lines respectively. 



